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In closed quantum systems, strong randomness can localize many-body excitations, preventing 
ergodicity. An interesting consequence is that high energy excited states can exhibit quantum co¬ 
herent properties, such as symmetry protected topological (SPT) order, that otherwise only occur in 
equilibrium ground states. Here, we ask: which types of SPT orders can be realized in highly excited 
states of a many-body (MB) localized system? We argue that this question is equivalent to whether 
an SPT order can be realized in an exactly solvable lattice model of commuting projectors. This 
perspective enables a sharp definition of MB localizability. Using this criterion, it is straightforward 
to establish that whereas all bosonic SPTs in spatial dimensions d = 1, 2, 3 are MB localizable, 
chiral phases (e.g. quantum Hall fluids) are not. We also show that free fermion SPTs in d > 1 
(including topological insulators and superconductors) cannot be localized if interactions are weak. 
A key question is whether strong interactions can render them MB localizable, which we study in 
the context of a class of d = 2 topological superconductors. Using a decorated domain wall (DDW) 
approach we show that some phases in this class are MB localizable, when they correspond to 
bosonic SPT orders. However, a similar DDW approach faces a fatal obstruction to realizing certain 
intrinsically fermionic SPT orders, an issue we argue may persist beyond this specific construction. 


I. INTRODUCTION 

Topological phases often exhibit unconventional edge 
states with unusual properties - for example the chiral 
edges of the quantum Hall phases, Majorana end states of 
topological superconducting wires, and the helical edges 
of topological insulators. In equilibrium settings, such 
topological edge states are protected by a bulk gap at 
zero temperature, but are, strictly speaking, destroyed at 
any finite temperature due to thermally activated excita¬ 
tions. An interesting exception arises in well-isolated and 
strongly disordered systems in which all excitations are 
localized, preventing the system from reaching thermal 
equilibrium - a phenomena known as many-body local¬ 
ization (MBL) PjEl This raises the intriguing possibility of 
harnessing coherent topological edge states in arbitrarily 
“hot” systems without the need for careful coolingPHU 

The possibility of MBL protection of topological edge 
states was recently explored in certain one-dimensional 
topological superconducting wires 7 and symmetry pro¬ 
tected phases^, as well as in the two dimensional 
symmetry-protected phase with Ising-like Z 2 symmetrj^ 
(Levin-Gu model—). Additionally, MBL was shown to 
protect intrinsic topological order in toric code type mod¬ 
els of discrete gauge theories.^® These studies have shown 
that the excited states of MBL systems have the prop¬ 
erties of quantum ground-states, including low entangle¬ 
ment and certain quantum orders. In this paper we would 
like to ask the converse question: which properties of a 
quantum ground-state be extended to all excited states 
in an MBL system? In particular, we aim to understand 
the scope of possible topological distinctions among dif¬ 
ferent MBL phases, and whether ground-state topology 
may be protected from thermal decoherence by MBL. 

To this end, we develop a more systematic understand¬ 
ing of the localizablity of topological phases, focusing 
primarily on “integer” topological- and symmetry pro¬ 


tected topological phases (SPTs). Building on ideas of 
Refs. (_5. Til and m we use the existence of a complete 
set of local conserved quantities as the defining feature 
of an MBL phase, and construct a sharp criterion for 
whether a phase is many-body localizable (MB localiz¬ 
able). Specifically, we define a phase to be MB localizable 
if it can be continuously connected (i.e. without encoun¬ 
tering a delocalization transition, or equivalently by a 
finite depth unitary circuit) to one whose entire eigen- 
spectrum is completely characterized by a set of strictly 
local conserved quantities. Here, we consider only “full” 
MBL in which the entire spectrum belongs to the same 
MBL phase. In practice, this criterion for MB localiz¬ 
ability is equivalent to asking whether a phase may be 
represented by an exactly solvable, local commuting pro¬ 
jector Hamiltonian. 

Apart from the connection to MBL, there are several 
motivations for constructing such exactly solvable lat¬ 
tice models for topological phases. Firstly, exactly solv¬ 
able lattice models play an important foundational role 
in the conceptual development of topol ogica l and symme¬ 
try protected topological (SPT) phased® serving as 
a proof-of-principle for the existence of new topological 
phases, and elucidating their properties through simple 
exact calculations. Moreover, the question of MB local¬ 
izability is tightly connected to questions about which 
phases of matter can be efficiently represented by a clas¬ 
sical computer. Namely, all eigenstates of a system de¬ 
scribed by a local commuting projector Hamiltonian can 
be represented as tensor product states (TPS’s) with fi¬ 
nite bond dimension! 1 ^ Conversely, given a TPS descrip¬ 
tion one can essentially always construct an exactly solv¬ 
able parent Hamiltonian (given a certain “injectivity” 
conditioipUBH). Lastly, the construction of exactly solv¬ 
able models may assist in the efficient simulation of in¬ 
teracting SPT edge states. Typically, simulating the d 
dimensional edge of a d + 1 dimensional SPT requires 
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simulating the full d + 1 dimensional system, due to 
the anomalous properties of the edge, adding substantial 
computational overhead and essentially prohibiting the 
direct simulation of surface states of 3D SPTs. However, 
the construction of solvable models for bosonic SPTs re¬ 
vealed the possibility of representing the d dimensional 
edge of a d + 1 dimensional bosonic SPT as a purely 
d-dimensional lattice model, just with a slightly uncon¬ 
ventional symmetry action. If possible, analogous con¬ 
structions for fermionic SPTs would be highly desirable, 
permitting for example the efficient simulation of inter¬ 
acting quantum spin Hall edges or 3D topological insu¬ 
lator surface states. 

With these motivations in mind and a sharp crite¬ 
rion for MB localizability in hand, we show that chiral 
phases such as integer and fractional quantum Hall sys¬ 
tems are fundamentally not MB localizable. We then 
go on to investigate the MB localizability of non-chiral 
higher dimensional SPT phases. By suitably modi fying 
the cohomology 23 24 and Walker-Wang constructions^ 8 25 
we show that all bosonic SPTs are localizable in all phys¬ 
ical dimensions, d = 1,2,3. 

Next, we turn to the localization of fermionic phases. 
Unlike bosonic SPT phases, fermionic SPTs permit a 
noninteracting description in terms of topological band 
theory. However, we establish a fundamental obstruc¬ 
tion to localizing such topological bands while preserv¬ 
ing symmetry; this generalizes the well known Wannier 
obstruction for quantum Hall phaseJ^Sl and the related 
time-reversal symmetry-based obstruction for topologi¬ 
cal insulators!^ We then argue, that the same obstruc¬ 
tion is not necessarily present in interacting fermion sys¬ 
tems - and construct a particular example of a local¬ 
izable interacting SPT of fermions using the decorated 
domain wall (DDW) approach of Ref. [23 Interestingly, 
this DDW construction fails for most fermionic SPTs ex¬ 
cept for a special set that are topologically equivalent to 
bosonic SPTs (once interactions are present). A related 
obstruction is present for Walker-Wang model based con¬ 
structions of fermionic SPTs - suggesting the interesting 
possibility that certain intrinsically fermionic SPTs fun¬ 
damentally resist localization in dimensions higher than 
one. 


II. DEFINING MANY-BODY LOCALIZABILITY 

Our first step is to formalize the notion of MB local¬ 
izability. A defining feature of MBL systems is local 
integrability: 5 1112 A system is defined as locally inte¬ 
grate if (almost) all its eigenstates can be labeled by an 
extensive set of well localized conserved quantities or “in¬ 
tegrals of motion”. Here, well-localized near a position r 
implies that the overlap between the integrals of motion 
and microscopic lattice degrees of freedom decays expo¬ 
nentially with distance from r. The definition also allows 
for a measure zero (in the limit of infinite system size) 
set of “resonances”, states which are not eigenstates of 


these integrals of motion. 5 -^ 

Familiar equilibrium phases of matter such as sym¬ 
metry breaking phases (ferromagnets, superfluids, etc...) 
are typically defined as the equivalence classes of quan¬ 
tum ground-states that can be continuously deformed 
into a representative “RG fixed point” product state 
without passing through a phase transition. By anal¬ 
ogy, it is natural to define an MBL phase as one whose 
Hamiltonian can be deformed into a strictly locally inte- 
grable Hamiltonian (one with a complete set of conserved 
quantities that have strictly bounded extent with no ex¬ 
ponential tails) without passing through a delocalization 
transition. 5 Here, we make the natural (though not rigor¬ 
ously proven) assumption that any system with exponen¬ 
tially well localized conserved quantities can be smoothly 
connected (e.g. to arbitrary precision by a finite depth 
unitary circuit) to a strictly localized one. 

As a concrete example, consider a Heisenberg chain in 
strong random field: L/xxx = — (hi&i + ■ <Ti+i). 

Here, for var (h) J 2 , the system is in an MBL 

phase with eigenstates labeled by emergent integrals of 
motiorpSl, rf « erf + ..., where (...) indicates small 
“dressing” of erf by other nearby spin operators (with 
coefficients that decay exponentially with distance from 
site i). It is quite natural to expect these dressed inte¬ 
grals of motion, r? can be mapped into the bare physical 
spins erf to exponential-in-depth accuracy using a finite 
depth unitary circuit. Hence, we can define the MBL 
paramagnetic phase of this model by those states that 
are connected (e.g. to desired accuracy by a finite depth 
unitary circuit) to the strictly local product eigenstates 
of Hxxx with J = 0, in exact analogy to how ordinary 
clean, T = 0 paramagnetic or symmetry broken phase 
are defined by adiabatic connection to a representative 
unentangled product state. 

Then, MBL phases may be represented as equivalence 
classes of models that are unitarily equivalent to a repre¬ 
sentative set of strictly local conserved quantities. With¬ 
out loss of generality we may take these conserved quan¬ 
tities as eigenvalues of a complete set of commuting pro¬ 
jector operators, {n r a }, that have bounded support in 
the region of position r, and which satisfy n 2 a = n ra . 
From this set of of projectors, we may construct a family 
of “fixed point” Hamiltonians for the MBL phase: 

H[X] = Y. X ^ U r,a (!) 

r,a=l...D 

Here a = 1... D is a band-index, and to achieve a sta¬ 
ble MBL phase {A r 0 } should be random set of strongly 
random coefficients chosen such that the variance in of 
A r ,a for nearby r is much larger than their typical value. 
The latter property ensures that the many-body local¬ 
ized character of the system is perturbatively stable to 
weak breaking of the strictly local integrability. We de¬ 
fine any system that is continuously connected to such 
an fixed point MBL Hamiltonian as many-body (MB) 
localizable. 

Then, a systematic understanding of MBL phases can 








3 


be obtained from studying the simpler fixed point Hamil¬ 
tonians H[ A]. In this paper, we are interested in the case 
where each eigenstate, |{n r , a }), of H[ A] has the prop¬ 
erties of a topological or symmetry protected topologi¬ 
cal ground-state, where n r a |{n r a }) = n r ,a|{n r ,a}) and 
n r ,a G {0,1}. For simplicity, we restrict our attention to 
“integer” topological states without fractionalized excita¬ 
tions. These could include “intrinsic” topological phases, 
or symmetry protected topological (SPT) phases. In the 
latter case, we further demand that the projectors be 
invariant under the action of the symmetry group G pro¬ 
tecting the topological phase. 

Though antiunitary (e.g. time-reversal) and continu¬ 
ous (e.g. charge-conservation) symmetries are more com¬ 
mon for electronic systems, these present additional the¬ 
oretical complications, and for this reason, unless other¬ 
wise specified, we will limit our analysis to SPTs pro¬ 
tected by discrete, unitary symmetry groups. 

Note, also, that our definition allows only for “full” 
many-body localization, in which the entire many-body 
spectrum is labeled by exponentially well localized inte¬ 
grals of motion. This leaves aside more exotic scenarios, 
such as the possibility of SPT orders that can only arise 
systems with a a partially MB localized spectrum, sepa¬ 
rated by thermal eigenstates by a MB mobility edge. 

The above arguments show that the question of which 
topological phases may be protected by MBL is equiva¬ 
lent to the question of which topological phases may be 
represented by commuting projector Hamiltonians. In¬ 
terestingly, these issues are closely connected to quantum 
information based questions about how to efficiently rep¬ 
resent quantum systems on a classical computer. Specif¬ 
ically, since the ground-state (in fact any eigenstate 30 ) 
of a commuting projector Hamiltonian can be repre¬ 
sented by a tensor product state (TPS) wave-function 
with bounded bond dimension, any MB localizable sys¬ 
tem permits an efficient TPS representation. The con¬ 
verse is also true so long as the tensors satisfy a certain 
injectivity conditioners 


III. PRELIMINARIES 

We are now in a position to understand whether or 
not a large class of topological phases permit MB local¬ 
ization. We begin with some simple examples that are 
relatively straightforward extensions of previously devel¬ 
oped methods, and then turn to the less well developed 
question of whether fermionic SPT phases may be MB 
localized. 


A. Chiral Phases 

In two dimensions there are a variety of chiral phases 
that are intrinsically topological (in the sense of not re¬ 
quiring symmetry protection) but not fractionalized (i.e. 


which have unique ground-state on a torus) - includ¬ 
ing integer quantum Hall systems and chiral topologi¬ 
cal superconductors. However, as pointed out in Ref. [151 
and 133 there is a fundamental obstruction to obtain¬ 
ing a commuting projector Hamiltonian of the form in 
Eq.|I]for topological phases with chiral edge modes. Chi¬ 
ral phases are partially characterized by a chiral central 
charge c c that determines the imbalance between left and 
right moving edge modes. Applying a temperature differ¬ 
ence AT across a finite sample of chiral phase produces 
a thermal Hall current density jn = c c yf yf TAT. While 
only the edge carries uniform heat current, this must be 
accomplished by having circulating thermal currents in 
the bulk, which cancel to give no uniform contribution 
except at the edge. However, the heat current density 
j_E(r, n; r' , m ) ~ [n r / im , n rj „] = 0 vanishes identically for 
commuting projector Hamiltonians. In the present con¬ 
text, this proof demonstrates that chiral phases cannot 
be MB localizedP 0j 3H Furthermore, recent worlPstrongly 
suggests that the chiral obstruction to localization also 
inevitably leads to ergodicity and thermalization in in¬ 
teracting chiral phases. 

While chiral edge modes present a fundamental ob¬ 
stacle to MB localization, one might expect that SPT 
phases, which by definition are smoothly connected (in 
the absence of symmetry) to trivially localized insual- 
tors and hence have non-chiral edges, may be MB local¬ 
ized. Indeed, we will see that many (though likely not 
all) SPTs are MB localizable. 


B. Bosonic SPTs 

All known SPT phases in d = 1,2,3, arising from 
bosonic degrees of freedom are MB localizable. Specifi¬ 
cally, group cohomology methods used to classify bosonic 
SPTs with discrete, unitary symmetries enable an ex¬ 
plicit construction of commuting projector Hamiltonians 
with bosonic SPT grounds statesP^EH As written, these 
commuting projectors form only an incomplete set, which 
uniquely specify an SPT ground-state but leave a highly 
degenerate excited state spectrum. However, as shown 
in Appendix |A 1[ these cohomology models can be suit¬ 
ably extended to include a complete set of projectors 
to produce a stable MBL system for which the entire 
spectrum of eigenstates take the form of SPT ground- 
states. The cohomology classification exhausts all known 
bosonic SPT phases in 2D, but misses a 3D bosonic SPT 
protected by time-reversal.^ This beyond cohomology 
state can be characterized by a symmetry respecting sur¬ 
face topological order (STO) that takes the form of a 
twisted Z 2 gauge theory in which all fractional parti¬ 
cles are fermions. This STO forms the basis of a differ¬ 
ent set of commutin g pro jector lattice models - Walker- 
Wang (WW) models^^Sl _ w hose uniform ground-state 
3D bosonic SPT phases with a given STO. Together, 
these constructions demonstrate that any bosonic SPT 
in arbitrary dimension is MB localizable. 
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Type of Phase 

Spatial dimension (d) 

MBL? 

Model type 

Chiral Phases 

2D 

N 


Boson SPTs 

ID,2D 

3D 

Y 

Y 

cohomology 

cohomology or Walker-Wang 

Free Fermion SPT 

ID 

2D,3D 

Y 

N 

noninteracting Kitaev chains 

Interacting Fermion SPTs 

2D 

Mixed? 

?? 


TABLE I. MB localizability of various integer topological phases - (‘Type of phase’ column) listed according to number 
of spatial dimensions, d , and whether or not the phase is MB localizable (‘MBL?’ column - ‘Y/N’ indicates Yes/No). If the 
phase is MB localizable - an example of the type of commuting projector model is given (‘Model type’ column). 


C. Free fermion SPTs 

Since noninteracting bosons simply condense to form 
a gapless superfluid, bosonic SPT phases require strong 
interactions to obtain an insulating state. By contrast, a 
non-interacting fermions can achieve an insulating SPT 
state simply by filling topological bands. Local commut¬ 
ing projector Hamiltonians of one-dimensional fermion 
SPTs are possible purely within such a free fermion de¬ 
scription, as fam ously exemplified by Kitaev’s supercon¬ 
ducting wire!- 1 ®! Generalization of Kitaev’s construction 
may be used to establish the MB localizability of all 
known one dimensional SPTs with various symmetries. 

However, in dimensions larger than one, free fermions 
filling SPT bands turns out to be a poor starting point for 
MBL, due to a fundamental obstruction to fully localizing 
the orbitals of topological bands. Such obstructions are 
well kno wn for chiral and time-reversal protected SPT 
phases,^ 2 ^, and also can be shown to extend (see be¬ 
low) to general fermion SPT bands with discrete unitary 
symmetry. This obstruction renders topological bands 
perturbatively unstable to thermalization upon the in¬ 
clusion of interactions.® 

To establish this free fermion obstruction for SPTs 
with unitary symmetries, consider a finite two- or three- 
dimensional system of noninteracting fermions with open 
boundaries described by commuting projector Hamilto¬ 
nian, ([!), with onsite unitary symmetry group G. Fur¬ 
ther, take uniform couplings, A r , a = —A for occupied 
bands (a = 1... n occ ) and zero for the remaining bands 
(a = n occ + 1... N). A finite system with open bound¬ 
aries can be chosen to have highly degenerate set of dis¬ 
persionless edge states by omitting projectors that in¬ 
tersect the boundary. We can remove this degeneracy 
(except perhaps a finite residual degeneracy required by 
symmetry) by adding noninteracting symmetry preserv¬ 
ing edge perturbations of characteristic energy scale V. 
For f«A, the edge bands reside inside the bulk gap. 
By contrast, for a nontrivial fermion SPT, the edge bands 
of a nontrivial SPT must disperse across the bulk gap and 
connect to gapped bulk modes in order to inherit the 
anomalous symmetry properties required for a nontrivial 
SPT edge state (see Appendix [b] for a formal proof). 
Hence, we see that free fermion commuting projector 


Hamiltonian can only describe trivial (not SPT) phases 
in dimensions higher than one. 

While these simple arguments rule out the MB lo¬ 
calization of higher dimensional non-interacting fermion 
SPTs, they leave open the question of whether interacting 
fermion SPTs can be MB localized, to which we now turn. 
Adding weak interactions to a non-interacting fermion 
“band” SPT is knowrPto further destabilize MB localiza¬ 
tion and lead to ergodicity and thermalization. Despite 
this, we will encounter an example in which strong, non- 
perturbative interaction effects can instead enable MB lo¬ 
calization. More generically, we will also give arguments 
that certain phases cannot be MB localized, even with 
arbitrary interactions. These results are summarized in 
Table U 

IV. 2D INTERACTING FERMION SPTS 

The key feature that distinguishes fermion systems 
from purely bosonic ones, is that any local operator con¬ 
tains an even number of fermions - implying that fermion 
parity is always conserved. Hence, given a symmetry 
group G, for a fermion system we should really consider 
a larger “symmetry” group G x Zf (assuming that Z^ 
is not already contained as a subroup of G) - where Zf 
is the Z 2 -group generated by the fermion parity operator 
Pp = (—1)™ F and hf is the fermion number operator. 
This special property of fermionic systems will turn out 
to have important implications to their MB localizability. 

In two dimensions, it is then always possible to intro¬ 
duce a 7r-flux defect such that fermions obtain a phase of 
(—1) upon encircling the 7r-flux. In a topological super¬ 
conductor, this 7r-flux is simply a superconducting vor¬ 
tex. For 2D fermion SPTs (f-SPTs), the symmetry prop¬ 
erties of the 7r-flux provide an important characterization 
of the phase, and delineate three distinct types of SPT 
phases: 

1. Projective f-SPTs - in which the 7r-flux con¬ 
tains symmetry-protected degeneracy associated 
with fermionic zero-modes on which symmetry acts 
projectively, and 

2. Fractional Charge f-SPTs - in which the 7r-flux 
is non-degenerate but carries fractional symmetry 
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charge. 


3. Bosonizable f-SPTs - in which the 7r-flux is in¬ 
variant under symmetry, 

Starting with a gapless symmetry preserving edge of a 
bosonizable f-SPT, it is possible to produce an energy 
gap for all fermion excitations at the edge without break¬ 
ing symmetry by condensing the 7r-flux at the edge (i.e. 
allow fluxes to freely tunnel in and out of the edge). This 
procedure involves strongly fluctuating phase-degrees of 
freedom - and hence requires interactions. Hence, an in¬ 
teracting bosonizable f-SPTs is topologically equivalent 
to a purely bosonic SPT accompanied by “spectator” 
fermions that form a trivial (Mott) insulator and do not 
play an essential role in producing the SPT order. An¬ 
other way to see this is to imagine gauging fermion parity 
(i.e. coupling fermions to a dynamical Z 2 gauge field). In 
the gauged theory, the 7r-flux is a gapped bulk quasipar¬ 
ticle, often called a vison, rather than a confined defect 
(that must be induced by an external field) as in the un¬ 
gauged theory. In a gauged bosonizable SPT phase, the 
vison is bosonic and has trivial symmetry properties - 
and can be condensed to confine the Z 2 gauge charges 
(i.e. the fermions) without disrupting the SPT order. 
The resulting phase is a purely bosonic SPT phase with 
no fermion excitations, but with the same SPT properties 
as the ungauged fermion SPT. 

In contrast, fermions play an essential role in establish¬ 
ing the SPT order when the 7r-flux has nontrivial sym¬ 
metry properties. For projective f-SPTs, the 7r-flux has 
symmetry protected degeneracy associated with bound 
zero-energy fermion modes. Consequently - condensing 
the 7r-flux at the edge (or gauging fermion parity and con¬ 
densing the vison) proliferates these gapless modes and 
does not produce an energy gap for fermion excitations. 
Similarly, in fractional charge SPTs the 7r-flux cannot 
be condensed without breaking symmetry, since the frac¬ 
tional charge of the 7r-flux cannot be screened away by 
the integer charged fermions. 

We will be able to show that all bosonizable f-SPTs are 
MB localizable with interactions - despite having a free 
fermion obstruction that prevents their localization in 
the absence of interactions. In contrast, projective f-SPT 
order presents an obstacle for MB localization, which we 
argue may potentially be fundamental, and present for 
arbitrary interactions. At this time, we are not able to 
make a definitive statement regarding fractional charge f- 
SPTs. Though exactly solvable models for systems with 
fractional charge f-SPT ground states have been obtained 
from the supercohomology construction^, these models 
operate on a constrained Hilbert space and it is presently 
unclear whether this construction yields a complete set of 
projectors that label all excited states. We leave this as 
a subject for future work (see also the discussion section 
below). 


V. OBSTACLE TO MB LOCALIZING CERTAIN 
INTERACTING 2D FERMION SPTS 

Let us give an example of how a non-chiral topological 
phase of fermions presents an obstruction to being de¬ 
scribed by a Hamiltonian composed of local commuting 
projectors. We will then extract general lessons for which 
phases may suffer the same obstruction. 

We consider a simple topological phase of 2D inter¬ 
acting fermions which can be viewed as a condensate 
of decorated domain walls (as described below) and use 
this picture to try to construct the commuting projector 
Hamiltonian. This approach will be shown to fail - which 
then raises the question of whether this is a fundamen¬ 
tal obstruction, or simply the limitation of the approach 
adopted. We then argue that it is indeed the former, by 
first showing that a commuting projector Hamiltonian for 
such a phase necessarily implies that it is a condensate 
of decorated domain walls, and any such condensate, in 
the context of the fermion SPT we are considering, will 
necessarily produce correlations that decay at least ex¬ 
ponentially. This poses a contradiction since correlation 
functions of separated local operators must have vanish 
in the zero correlation length ground states of commuting 
projector Hamiltonians. 

Consider a 2D topological superconductor with a 
unitary Z 2 symmetry (in addition to the ever-present 
fermion parity conservation “symmetry”), which has a 
classification labeled by integer topological invariant A f. 
The J\f = 1 SPT phase is obtained by combining oppo¬ 
site chiral superconductors, one a p x + ip y and the other 
a p x — ip y superconductors where the Z 2 symmetry acts 
trivially on the fermions in the p x +ip y layers, and gives a 
phase (—1) to those in th ep x —ip y layer. The edge consist 
of a right and left moving chiral Majorana fermions "f s 
where s = R, L labels right and left movers respectively, 
described by the continuum Hamiltonian: 

-Hedge = -iv J dx'yja^ s ,d x 'y s > (2) 

the mass term H m = if dxM^jL changes sign under 
the Z 2 symmetry and is hence forbidden. As for any Z 2 
symmetric state, this M = 1 state can be viewed as a 
string-condensate of Z 2 domain walls in the symmetry¬ 
breaking mass, M. However, unlike a trivial paramag- 
net, the domain walls of the Z 2 phase are topologically 
non-trivial, and can be thought of as being “decorated” 
with ID TSC chains that possess unpaired Majorana end 
stated. This decoration is revealed by performing a 
gedanken experiment in which the Z 2 symmetry is bro¬ 
ken in two opposite ways in the top and bottom halves of 
the system. The resulting domain wall terminates at the 
boundary where it is associated with a change in sign of 
the mass M(x) = Mq ■ sign(a:). As can be readily verified 
from Eq. [2] such a domain wall carries binds a Majorana 
zero mode: 7 o = / dx e ^ [7 l{x) + 7 ^( 0 :)], near x = 0 , 

where £ « 
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Hence we can view the 2D Af = 1 TSC state as a fluctu¬ 
ating paramagnetic state of Z 2 “spins” whose DWs are 
“decorated” with ID Majorana chains. Since the bulk 
domains have closed boundaries, and the decorating ID 
SPT has an energy gap, so too does the bulk of the deco¬ 
rated domain wall (DDW) state. Related DDW construc¬ 
tions have been discussed for bosonic SPT phase^D. Here 
we attempt an analogous construction for the fermionic 
2D TSC, to generate parent Hamiltonians of local com¬ 
muting projectors. We will find that this fermionic DDW 
construction suffers a fatal flaw. 


A. Explicit DDW construction attempt 

To implement the DDW construction, we introduce 
bosonic degrees of freedom that transform as doublets 
under the Z 2 symmetry. Specifically, we consider Ising 
spins Tp = ±1 that reside on plaquettes, P, of a honey¬ 
comb, and complex fermions with annihilation operators 
c r on the sites r of the honeycomb. The Z 2 subgroup of 
the symmetry flips the Ising spins but does not effect the 
fermionic degrees of freedom. 

A protypical Z 2 symmetric state can be expressed as a 
superposition over all Ising spin-configurations, or equiv¬ 
alently as a equal superposition of domain wall (DW) 
configurations C that divide t z = ±1 domains: 

|* PM ) = ^ $>[£],« (3) 

where {^[C])} are DW-configuration dependent wave- 
functions for the fermionic degrees of freedom. We can 
produce a nontrivial 2D SPT phase by choosing | ^[C]) 
to a state in which fermions residing on a Z 2 domain 
wall form a ID TSC chain, while fermions residing inside 
ordered Z 2 domains form a trivial insulating state. 

This “decorated domain wall” (DDW) approach has 
previously been used to obtain (non-commuting projec¬ 
tor) lattice models of bosonic SPTsP 8 Here, we general¬ 
ize this construction by decorating Ising DWs with ID 
fermionic SPTs, and by modifying the terms that prolif¬ 
erate DWs to obtain a commuting projector model. 

The Hamiltonian consists of projectors that decorate 
DWs of Tp with Kitaev chains, and transverse field terms 
that flip Tp and rearrange the fermions to preserve the 
DW decorations. The domain wall decoration is accom¬ 
plished by: 

H decoiate = J2 U cH F {C) (4) 

c 

where indicates a sum over all r configurations C, 
and H F (C) takes the form of a zero correlation length 
Hamiltonian of a Kitaev TSC chain for fermion links re¬ 
siding on the domain walls of r and is the Hamiltonian of 
a trivial product state insulator for sites within ordered 
domains of r. To explicitly write the decoration terms, it 



FIG. 1. Schematic representation of fermionic DDW 
lattice model - The model consists of Ising spins, r sit 
on green sites, which form triangular lattice, and complex 
fermions, c r at sites, r of a honeycomb, which may be de¬ 
composed into two real Majorana fermions \ ( re d dot) and 77 
(blue dot). 



FIG. 2. Domain wall condensation (“transverse field”) 
terms - (a) Projector with no DWs. (b,c) Schematics of 

selected projectors that attach Kitaev chains to DWs. Purple 
link with arrow from to 77 r > indicates a iXrilr 1 term for 
fermions. 

is convenient to represent the complex fermions by Ma¬ 
jorana fields Xr = c r + cl and r] r = — i ( c r — cl) , such 
that c r = (Xr + if]r)/ 2, then: 

H f {C) = - ^2 iiir+dXr ~ iXrT]r ( 5 ) 

r£DW rj!DW 

Here, one needs to choose an (arbitrary) orientation con¬ 
vention for the DWs, for example the DWs can be ori¬ 
ented clockwise around the r 2 = +1 domains. 

To restore the Z 2 symmetry, one needs to proliferate 
DWs in the r spins, while preserving the DW decoration. 
This is accomplished by the generalized “transverse-field” 
terms: 

-Htf = T,J2 T P U c\F diP) )(F c \ (6) 

p c 

Here, )T) C indicates a sum over all r configurations C, 

| Fc) is the fermionic wave-function corresponding to the 
ground-state of Eq. [4] with r’s in configuration C, and 
C (P) is the spin configuration obtained from C by apply¬ 
ing Tp to the spin on plaquette P. Importantly, despite 
superficially appearing to depend on the global spin con¬ 
figuration C, for fixed P, the generalized transverse-field 
terms in Eq. [6] actually involve only a bounded set of op¬ 
erators in the vicinity of P. This is because the ground 
state of the decoration terms in Eq. [4] are determined by 
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purely local rules, and hence \Fq^)(Fc(p)\ involves only 
those fermionic degrees on the sites bordering plaquette 
P. 


1. Fermion Parity Obstruction 

While the above construction appears to yield an ex¬ 
actly solvable lattice model of a 2D interacting TSC, 
it actually contains a fatal flaw. Specifically, a com¬ 
plication arises due to the fact that the ground-state 
of a single Kitaev superconducting chain with periodic 
boundary conditions has odd fermion parity compared 
to the trivial ground state. To see this, consider the 
fermion parity operator along a closed loop of sites F: 
Pf = IXrer iXrVr = — Ilrer iVrXr+i- In the trivial 
ground state, iXrVr = 1 for each site, and in the ID TSC 
ground state with periodic boundary conditions along F, 
ir] r Xr+i = 1, demonstrating that Pp is opposite for the 
trivial and topological loops. This seemingly innocuous 
observation will turn out to present a serious obstacle to 
MB localizing the AT = 1 phase. 

Clearly, the parity switching feature presents a prob¬ 
lem for the DDW construction with J\f = 1 since, here, 
the desired TF operators \Fc(p))(Pc\ actually change 

fermion parity whenever C(P) and C differ by an odd 
number of domain walls. Such fermionic operators are 
highly nonlocal (they anticommute with other fermion 
operators even arbitrarily far from plaquette P) and 
hence are not permitted in a local Hamiltonian. At this 
point we may conclude that perhaps the DDW approach 
is not inadequate in this particular situation or may need 
to be modified e.g. by considering condensates of domain 
walls that enclose n flux, and hence even fermion parity. 
We show below that neither of these is a viable option. 
First we argue that a condensate of n flux domain walls 
has correlations that decay at least exponentially with 
distance. Such correlations should be absent in ground 
states of local projectors. On the other hand, we argue 
that such a condensate is mandated by the assumption 
of a parent Hamiltonian of local commuting projectors, 
which forces the desired contradiction. 


2. Finite Correlations in the Decorated Domain Wall State 

The difficulty of nontrivial fermion parity of the Ma- 
jorana chains decorating the domain walls can be cir¬ 
cumvented by twisting the boundary conditions to add 
a 7r-flux Eq. [6j but we will see that restoring locality 
to the DDW Hamiltonian in this way comes at the ex¬ 
pense of spoiling the exact solvability of the model. To 
see this, consider the operator (D(r) acting on a minimal 
plaquette at position r that adds a flux-free domain with 
appropriately decorated domain wall. Specifically, acting 
on |^>[C]), 0( r) either creates a fluxless loop if the pla¬ 
quette containing r resides inside a connected domain of 


the configuration C, or grows an existing domain (with¬ 
out adding additional flux), if the plaquette borders a 
domain of C (see Fig. [3]). 

The expectation value of two such operators 
(O(r)O(r')) does not factorize into (0(r))(0(r')) even 
for widely separated r, r'. To see this, consider the 
loop configuration, C, shown in Fig. [3] Acting on this 
configuration O(r) merges two domains (each with n- 
flux) without adding any flux - producing a domain with 
2-7T ~ 0 flux. Such fluxless domains are not present in the 
ground-state |^ 2 dspt), implying ('F 2 dspt|C , ( i ')I 1 / ; [C]) = 
0. On the other hand, acting on | if>[C}) with both 0(r) 
and O(r') merges three domains together - and hence 
preserves the odd flux through the domains, implying: 
( v I / 2DSPT|C>(i')C , (r , )l' i /’P]) i 1 0. 

In general, 

AT/ 

C{r,r') = (O(r)O(r')) - (O(r))<0(r')> « ^ 

(7) 

where N' c is the number of configurations of the sort 
shown in Fig. [3] These special configurations require 
r and r' to be connected by a single symmetry domain. 
Since the SPT ground-state (by definition) must preserve 
symmetry. Then, the ratio is proportional to the 
probability that a single domain percolates from r to r'. 
In any gapped symmetry preserving phase, this proba- 
biltiy decays exponentially with separation |r — r'| with 
a non-zero characteristic lengthscale £. This non-zero 
correlation length, £, contradicts the assumption that the 
phase can be described by a strictly local commuting pro¬ 
jector Hamiltonian whose correlation functions factorize 
exactly beyond some finite lengthscale. We note that al¬ 
ternative methods for circumventing this fermion parity 
obstruction, e.g. by adding extra fermionic degrees of 
freedom to the DW loops, are topologically equivalent to 
condensing 7r-flux loops and should suffer from the same 
non-locality. 

These considerations strongly suggest that the require¬ 
ment that each symmetry domain contain a 7r-flux is in¬ 
compatible with the commuting projector form of the 
Hamiltonian, and that, consequently, there is a funda¬ 
mental obstruction to deforming this AT = 1 TSC to the 
strictly localized limit. This obstruction is present for 
any projective f-SPT where the ID SPT decorating the 
symmetry domain wall is topologically equivalent to a 
collection of Kitaev TSC chains. 

A potential shortcoming of the above argument is that 
we assume that it is possible to construct a flux-less in¬ 
sertion operator, O , that has non-vanishing two-point 
functions in some eigenstate. While there is no symme¬ 
try principle requiring {0{r)0{r')) to vanish in all eigen¬ 
states, it is conceivable that one may be able to construct 
a highly fine-tuned model for which this two-point func¬ 
tion vanishes in every eigenstate, possibly enabling a loop 
hole in our argument. 



c) 



a) 




n R 


Y>W)—<i—<i-(!> 

"Y—T—Y — Y—T 


FIG. 3. Non-zero correlation length - A schematic il¬ 
lustration of a DW configuration, C, for which the correlation 
function of the flux-less loop insertion operator O fail to fac¬ 
torize - implying non-zero correlation length and an obstacle 
to constructing a strictly local projector model. 

3. Necessity of the DDW perspective for J\f = 1 

At this point, one might suspect that the above prob¬ 
lems merely demonstrate a shortcoming of the DDW ap¬ 
proach for building solvable lattice models of fermionic 
SPT phases. However, we will argue that the DDW de¬ 
scription and its associated problems for creating a com¬ 
muting projector Hamiltonian are fundamental to the 
J\f = 1 state, and therefore unavoidable. 

We first show that, given a putative local commut¬ 
ing projector Hamiltonian for the Af = 1 TSC with Z 2 
symmetry, we can formally construct a family of DDW 
insertion operators that commute with the Hamiltonian. 
For the reasons outlined above, this then suggests a non¬ 
vanishing correlation length in contradiction with the as¬ 
sumption of a local commuting projector Hamiltonian. 

A DDW insertion operator acts with the Z 2 symmetry 
generator inside a closed connected region, A, creating an 
extra Z 2 domain, and also inserts a Kitaev chain along 
the boundary. To define a local action of symmetry con¬ 
sider acting with the generator of the Z 2 symmetry, g , 
in a bounded region, A , implemented by the operator: 
UA{g) = IlieA gi, where gi is the generator of group el¬ 
ement g acting on physical site i. Under the action of 
UA(g) the original projectors {n} map to a new set of 
commuting projectors: n ra = UA(g)W^ a U a(<?)^ (note 
that {IT} are also a commuting set of projectors). Pro¬ 
jectors residing in the interior of A , and in its complement 
A c are unchanged, but those whose support intersects the 
boundary of A, dA are generically altered. However, we 
may repair projectors overlapping dA by acting a string 
operator, WoA{g), with bounded support in the vicinity 
of dA , that maps the projectors transformed by £( 4 ( 3 ) 
back to the original projectors. Thus, for any local com¬ 
muting projector Hamiltonian, we may “pull back” the 


FIG. 4. Formal construction of DDW insertion oper¬ 
ator - (a) Given a commuting projector Hamiltonian, one 

can construct an operator that commutes with the Hamilto¬ 
nian and inserts a symmetry domain - i.e. which acts with 
a symmetry group generator, g, inside a region A (blue glob¬ 
ular region), and a string operator, Wqa, on its boundary, 
dA (hatched strip), (b) By examining the end of the opera¬ 
tor Wqa where it intersects a spatial boundary (orange star), 
and using the bulk-boundary correspondence, one can show 
that for projective fermionic SPTs, Wqa inserts a ID SPT. 
(c) Though independent, projector terms in the Hamiltonian 
(schematically indicated by green and orange squares) gener¬ 
ically overlap multiple physical sites (gray circles). 


local action of symmetry on a 2D domain, into a ID 
string operator on its boundary. 

This formal construction yeilds an infinite family of 
DW insertion operators, {££ 4 ( 3 ) = Wo a (g)UA ( 3 )}, 
which commute with the Hamiltonian. Hence, the 
ground-state (indeed all eigenstates) are eigenstates of 
{Da ( 3 )}, and must be regarded as “condensates” of these 
DW insertion operators. 

To complete the DDW interpretation, let us consider 
the end of string operator WoA{g), formally obtained by 
acting with ££4(3) on a semi-infinite system in the left 
half plane ({x < 0 ,y £ R}) with open boundary along 
x = 0, and open boundaries, and take A to be the lower 
half-plane {y < 0). By construction, the projectors that 
do not overlap the point (x = 0 ,y = 0), are left un¬ 
changed, the operator £>, 4 ( 3 ) only effects degrees of free¬ 
dom near x = 0 , y = 0 and thus creates a low-energy 
excitation on the boundary, i.e. ££ 1 ( 3 ) is the bound¬ 
ary operator that inserts a Z 2 symmetry twist (possibly 
dressed by some local boundary operators at the origin). 

Acting on a vortex (7r-flux) in the J\f = 1 2D TSC, the 
fermion parity operator anticommutes with the Z 2 sym¬ 
metry generator, i.e. PfDa^Pf = — ££ 1 ( 3 )- Since Pf 
commutes with the local action of Z 2 , this implies that 
PfW A(g)PF = —WA{g) - i-e. that the ID SPT deco¬ 
rating the “domain-wall” of the Z 2 symmetry develops 
odd fermion parity when encircling a 7r-flux. We have 
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already seen that a loop of ID TSC, has precisely this 
property, in fact, this property is sufficient to define the 
ID TSC phase - in the sense that any ID state with this 
property is topologically equivalent to that described by 
([ 5 ]). This shows that in any putative commuting projec¬ 
tor Hamiltonian, one can formally construct DDW in¬ 
sertion operators that commute with the Hamiltonian, 
i.e. that the DDW interpretation of the Af = 1 state is 
indispensable, and that the concomitant barriers to MB 
localization encountered in the previous section are po¬ 
tentially fundamental. Moreover, these arguments read¬ 
ily extend to other projective f-SPTs, which have similar 
DDW interpretations and face a common barrier to MB 
localization. 


B. DDW Construction for Af > 1 

While the DDW viewpoint is not strictly necessary for 
AT > 1, we may still attempt to construct solvable lat¬ 
tice models for these phases using the DDW method by 
extending the above model to include Af > 1 flavors of 
fermions, c rjQ , labeled by flavor index a = 1.. .Af. How¬ 
ever, as there is only one kind of non-trivial ID TSC wire, 
to implement the DDW construction for Af > 1, we must 
artificially extend the the symmetry group Z 2 —> Z 2 x G, 
where G is a (to be specified) internal, unitary symmetry. 
We can then decorate the Z 2 DWs by ID SPTs protected 
by symmetry group G. 

At first glance, the odd fermion parity problem for 
Af =1 does not affect the DDW models with even num¬ 
ber of fermion flavors Af', however, a related problem 
emerges for Af = 2. To see this note that, for Af = 2 
terms such as c ri ic rj 2 must be symmetry forbidden, oth¬ 
erwise it would be possible to deform the DW decora¬ 
tion into a trivial insulator without breaking symme¬ 
tries. Equivalently, the two different species of fermions, 
a = 1,2 must have distinct charges under the residual 
symmetry group, G, of the DW (recall that the DW con¬ 
figuration breaks the Z 2 subgroup of the total symme¬ 
try group Z 2 x G). Suppose we decorate the Ising DWs 
with ID fermion SPTs having periodic boundary con¬ 
ditions, as required to obtain a zero-correlation length 
model. Then the requirement that c rjQ transform differ¬ 
ently under G for a = 1, 2, implies that \Fc(p))(^C | have 

nontrivial G-charge whenever C(P) and C differ by an 
odd number of domain walls. 

For Af = 4, however, such obstructions need not occur. 
One suitable option is to take G = Z 2 x Z 2 , and give 
each of the four flavors a different pair of charges from 
the set from the set {(+, +), (+, —), (—,+), (—, —)} under 
the two Z 2 subgroups of G. For this choice, each of the 
four copies of Kitaev chains decorating the Ising domain 
walls has distinct transformation properties, ruling out 
the bilinear mass terms c r ^ a c rt b- However, the TF opera¬ 
tors \Fc( P j)(Fc\ all have trivial symmetry quantum num¬ 
bers as the product of charges of all four fermion species 
is trivial. Indeed, for this specific choice, one can read¬ 


ily verify that the DDW Hamiltonian, (4[6), produce the 
desired symmetry-preserving commuting projector form 
(see Appendix D 4 for details). 

Interestingly, the obstruction to localization is lifted 
precisely because the fermion phase becomes bosoniz- 
able at Af = 4. The underlying bosonic character of the 
phase obtained from the DDW construction with Af = 4 
can be seen by adding interaction terms n:=i Xr,a and 
n a= i Vr,a to the edge of the system. These gap out all 
fermionic excitations at the edge, but leave symmetry 
protected gapless bosonic modes (see Appendix |D 5| ). Al¬ 
ternatively, one can examine a 7r flux for the fermions 
and verify that it has trivial symmetry properties (see 
Appendix D5). Hence, we could remove the fermions 
entirely from the problem without effecting the SPT or¬ 
der by coupling the fermions to a Z 2 gauge field and 
condensing the gauge-flux to confine the fermionic ex¬ 
citations. In this construction, the choice of symmetry 
implementation that enables the DDW construction to 
work, also allows the interaction terms n;:= =1 Xr,a and 
Ha —1 Vr,a that can gap out the fermion degrees of free¬ 
dom at the edge - producing an essentially bosonic SPT. 
This construction then suggests that whether a fermion 
SPT is bosonizable and whether it is MB localizable are 
related issues. 


VI. GENERALIZATIONS AND DISCUSSION 

The above-described obstruction appears to rule out 
MB localization of 2D projective f-SPTs (DDW phases), 
and highlights a sharp distinction between bosonic SPT 
phases which can all be MB localized and permit exact 
TPS descriptions, and fermion SPTs for which can face 
obstructions to MB localization unless they are bosoniz¬ 
able. 

At this time, we are unable to decisively state whether 
the fractional charge f-SPT phases also resist MB local¬ 
ization. The simplest scenario is that all intrinsically 
fermionic (i.e. not bosonizable) SPT phases share a com¬ 
mon obstacle to MB localization, i.e. that the special 
properties of fermion parity “symmetry” make it impossi¬ 
ble to endow 7r-fluxes (or flux lines in 3D) with anomalous 
symmetry properties. Further suggestive evidence for 
this scenario comes from considering Walker-Wang type 
loop models for 3D fermion SPTs. Unlike their bosonic 
counterparts, which give an MB localizable family of pro¬ 
jector Hamiltonians, fermion Walker-Wang constructions 
actually produce a gauged version of the fermionic SPT 
order in which the fermions are emergent anyonic parti¬ 
cles coupled to a dynamical Z 2 gauge field. These loop 
models are only symmetry preserving in the zero-gauge- 
flux sector - since for a fermionic SPT the Z 2 gauge flux 
lines (7r-flux lines for the fermions) are necessarily gapless 
or break symmetry. Moreover, while it is possible “un¬ 
gauge” these models by coupling them to local fermion 
degrees of freedonPD, this ungauging action spoils the lo¬ 
cality of the terms in the Walker-Wang model and does 
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not produce a local commuting projector Hamiltonian. 
The failure of these Walker-Wang constructions due to 
the special properties of 7r-flux lines in fermion SPTs is 
consistent with a scenario in which fermion SPTs are in¬ 
trinsically not MB localizable. 

However, the supercohomology results of Ref. [331 sug¬ 
gests that the situation may be more complicated than a 
simple dichotomy between fermionic and bosonic phases, 
i.e. that only certain types of fermionic phases (e.g. pro¬ 
jective f-SPTs) face an obstruction to MB localization. 
In particular, exactly solvable models with a fermionic 
SPT ground states (specifically a fractional symmetry 
charge 2D f-SPTs) was constructed by supercohomology 
methods in Ref. [33]. While it is conceivable that this 
constructi on m ay be extended (e.g. along the lines of 
Appendix. A 11 to obtain an MB localizable Hamiltonian 
for all excited eigenstates, an explicit demonstration is 
absent. An added complication compared to the bosonic 
case is that the fermionic models operate in a restricted 
Hilbert space, which may present complications for con¬ 
structing a generically MB localizable model. These su¬ 
percohomology results suggest further subcategorization 
within fermion SPT phases - in which only projective 
f-SPTs suffer a barrier to MB localization. 


A third possible scenario is that only phases like the 
J\f = 1 2D topological superconductor phase, in which 
symmetry defects carry unpaired Majorana zero-modes, 
suffer an obstacle to localization. Distinguishing among 
these scenarios will be an important subject for future 
work. 


A second limitation of our current work is that we 
have so far considered only SPT phases protected by dis¬ 
crete, unitary symmetry groups. Our arguments against 
MB localization in projective f-SPTs can likely also be 
extended to antiunitary time-reversal symmetry, since, 
for local commuting projector Hamiltonians it is possi¬ 
ble to go to a tensor-product state description where one 
can define the local action of an antiunitary symmetry.^ 
Extending these ideas to SPT phases protected by con¬ 
tinuous symmetry group, however, appears problematic. 
For example, in 2D there is an integer quantum Hall 
phases of bos ons p rotected by a 17(1) charge conserva¬ 
tion symmetry)2EL 3 *JThe group cohomology approach also 
produces exactly solvable lattice Hamiltonians for these 
phases. However, these models utilize U(l) rotor de¬ 
grees of freedom with an unbounded Hilbert space on 
each site. The unbounded character of the on-site Hilbert 
space, inevitably allows for resonant transitions for any 
given pair of sites somewhere in the spectrum. A related 
conjecture, that continuous symmetry SPT phases resist 
localization, was posed for one-dimensional spin chains 
with continuous rotation symmetry®. Further support 
for this hypothesis for the case of continuous non-Abelian 
symmetries was obtained by more systematic analysis 3 -®. 
Whether the infinite on-site Hilbert space can be trun¬ 
cated without spoiling the exact solvability of the model, 
or alternatively whether the accidental resonances in the 
rotor models are sufficiently prevalent to completely spoil 


localization are subtle questions that we leave for future 
work. 

We remark, however, that, if possible, MB localiza¬ 
tion of these bosonic IQH states would have rather 
striking consequences: since the U(l) degree of freedom 
cannot be spontaneously broken at the ID boundary 
of a 2D system, due to Mermin-Wagner arguments, 
bulk localization of the bosonic integer quantum Hall 
phases could topologically guarantee the presence of 
gapless, ergodic boundary modes. This situation would 
be unlike that for discrete-symmetry SPT phases, 
whose symmetry can always be spontaneously broken 
at the edge, allowing a localized non-ergodic edge 
termination. These considerations hint at interesting 
connections between MB localizability and Mermin- 
Wagner type arguments that deserve further exploration. 
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SPT order. 


Appendix A: MB localizable models for bosonic 
SPTs 


1. Completing the cohomology models 

Exactly solvable lattice models with SPT ground- 
states, were constructed in Refs. E51 and f231 for all bosonic 
SPT phases classified by group cohomology. - i.e. those 
with discrete unitary on-site symmetry group G. The 
models involve an appropriate d-dimensional lattice (e.g. 
triangular for d = 2), with bosonic site-degrees of free¬ 
dom {|<?)i}, on which G acts as such that gi\g 2 ) = 
I <7i <72 ) - The construction starts from a trivial Hamilto¬ 
nian: H tliv = - 10*}(0,1 where |0») = E 3 eG Iff*)- 

The SPT Hamiltonian is then obtained as: 

J7o = -^n* = -^C7 i |0 i >(0 i |C// (Al) 


where [/* are unitary matrices diagonal in the 
| ... gigj ...) basis, with elements consisting of phases 
chosen from the group cohomology: 


n = E 

I g} 


riijfceA t'd+i (i) gii gj ) gk) 
ILjfcev Vd+i(i,gi,gj,gk) 


IMXMI 


(A2) 
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where Vd+i are phases obtained from a particular element 
(SPT phase) of the cohomology group 'H d+1 (G, 17(1)), 
which involve the g- values of nearest neighboring sites of 

i. 

While H 0 has a unique non-trivial SPT ground-state, 
the commuting projectors II^o are not complete for |G| > 
2 (where |G| is the number of group elements in G). This 
is because there are |G| — 1 states per site that are or¬ 
thogonal to Ui\0i). For |G| > 2, these |G| — 1 orthogonal 
states are all annihilated by H , implying that the excited 
states are highly degenerate. Hence, as written the model 
is not suitable for MBL. 

Fortunately, this issue can be remedied by extending 
n ii0 to include complete set of projectors such that every 
state in the spectrum behaves as an SPT ground-state 
in the same phase. Let us specialize to a finite Abelian 
group G, where we may conveniently work in the Fourier- 
transformed basis: |fci) = YhgeG e zk ^ Iff)■ For ex¬ 
ample G = 7jn there is a single generator g, and all group 
elements can be written as g n with n € 0 ... N — 1 , and 
I h) = y= J2 g eG e27likn \9 n ), and k <E {^}j=o,...jv- 1 - 

More generally, any finite Abelian group can be decom¬ 
posed as a product of a finite number of Z^v’s: G = 
'Ln 1 x ... • In this case k will be an m-component 



The Hamiltonians Hk = —^2iUi\ki){ki\U} are each 
related to H 0 by the unit-depth local unitary transfor¬ 
mation: Hk = M(k)HoM(k )^, with M(k) = JX,Mi(fc), 
and Mi{k) = J2 gi e lk{9i) \gi){gi\, implying that H k ’s 
ground-state is the same SPT phase as H 0 ’s for each 
k. To see this, note that Mj(fc) are diagonal in the 
| gi) basis and hence commute with Ui implying: H k = 
- Ei u i M i (&)|0*)(Oil-M/ {k)Uj =M(k)H 0 M(k) t. More¬ 
over, it is easy to check that this unitary transforma¬ 
tion preserves the commuting projector structure. Fi¬ 
nally, the projectors are manifestly invariant under the 
symmetry G, since \ki ) are eigenstates of the symmetry 
group generators whse eigenvalues have unit amplitude, 
and since the Ui are constructed from the cohomology to 
be G-invariant. 

Then, the commuting projectors form a complete 
set that span all excited states of the system, enabling one 
to construct MB localized family of fixed point Hamil¬ 
tonians, Eq. [l] for each SPT phase in the cohomology 
classification. 


Appendix B: Free Fermion Obstruction 

In the main text, we argued that the symmetry re¬ 
specting edge of an SPT must exhibit bands that disperse 
across the bulk gap, and that, since commuting projector 
Hamiltonians produce flat edge bands and since symme¬ 
try cannot be broken in a non-interacting fermion system 
it should be impossible to have a MB localizable SPT 
of noninteracting fermions. To formalize this intuition 


let us consider first the case of a 2D SPT with an on¬ 
site, unitary symmetry group G. This argument readily 
extends to 3D, by considering one dimensional cuts of 
the 2D boundary Brillouin zone. The most general non¬ 
interacting SPT edge contains equal number N, of left 
and right moving Majorana fermion modes 'Jr/l^ with 
a € { 1 ... N} (complex fermion zero modes are special 
case of N even). The left and right movers must cross at 
some point within the Brillouin zone, and without loss of 
generality we may assume that the crossings occur mo¬ 
mentum k = 0 along the edge. 

In order to have the edge bands lie within the bulk 
gap, the right movers and must convert into left movers 
at the Brillouin zone boundary k = ±— where a is the 
lattice spacing. This requires that at k = ir the left 
and right movers must intersect in groups that transform 
under the same representation of symmetry: 7 l/r,o, 

( Og) a h lL/R,b where O g is an orthogonal N x N repre¬ 
sentation of the group element g £ G. However, the 
physically consistent orthogonal representations of G are 
discrete (even for continuous symmetry groups like G(l)) 
and cannot continuously evolve as a function of edge mo¬ 
mentum k. Hence, we may choose a basis at k = 0 where 
left and right movers transform identically under symme¬ 
try, permitting the mass term: M + h.c. 

which trivially gaps the edge - implying that the bulk 
cannot be an SPT. 


Appendix C: Review of 2D Topological 
Superconductors with Symmetry 

The three classes of 2D fermion SPT phases (projec¬ 
tive, fractional charge, and bosonizable), described in the 
main text, are all neatly illustrated by a family of 2D 
topological superconductors (TSCs) with a unitary Z 2 
symmetry. These examples are closely related to those 
introduced in Ref. HH The Z 2 symmetry can act on, 
for example, orbital or layer degrees of freedom (while 
time-reversal symmetry is perhaps a more natural dis¬ 
crete symmetry in electron systems, its antiunitary char¬ 
acter unecessarily complicates the theoretical discussion). 
The SPT phases of these TSCs are labeled by an integer, 
J\f, and can be thought of as J\f copies of p + ip super¬ 
conductors made of spinless (or spin polarized) electrons 
that are invariant under Z 2 , and A f copies of p — ip super¬ 
conductors made of electrons that obtain a (— 1 ) phase 
under the Z 2 symmetry (i.e. have an integer unit of Z 2 
“charge”). 

The edge of the Af th phase contains A f right-moving 
Z 2 -neutral Majorana modes, 7 r i0i (a = 1.. .J\f) and J\f 
left-moving Z 2 -charged Majorana modes, 7 L,a- In non¬ 
interacting systems all backscattering terms * 7 _R, a 7 L,b are 
symmetry forbidden, and hence the gapless nature of the 
edge is ensured by symmetry for any integer A f. In in¬ 
teracting systems, however, phases whose M -value differs 
by a multiple of 8 become topological equivalent - indi¬ 
cating that the integer (Z) classification of interactions 
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collapses to a Z§ interacting classification.^HH 

To see how this works, let us consider the properties 
of a superconducting vortex (7r-flux) for various A f. 


1. Af = 1, Projective f-SPT 

Consider first, the minimal case Af = 1. A supercon¬ 
ducting vortex in a single p + ip superconductor carries 
an unpaired majorana zero-mode, y (and similarly for 
p — ip superconductors). Then, a vortex in the Af = 1 
TSC carries two zero inodes: a Z 2 -neutral mode 7 + from 
the p+ip layer and a Z 2 -charged mode 7 _ from the p— ip 
layer. This pair of zero modes defines a complex fermion 
/ = 7+ ~!^ 7 ~ and an associated two-level system with 
states |±) defined by: iy + y_|±) = ±|±), interrelated by 
|+) = /t|—). Since |+) = fi\ —}, the states |±) have op¬ 
posite fermion parity. On the other hand, the Z 2 symme¬ 
try generator interchanges these states: Z 2 : |±) —► |=f). 
To see this, denote the generator of Z 2 symmetry as g , 
then: g (ry+ 7-)5 = * 7 +(~ 7 -) (alternatively gfg = / t ). 

Moroever, when acting on the vortex zero modes, 
the Z 2 symmetry anticommutes with the ever-present 
fermion parity (Zf) “symmetry” generated by opera¬ 
tor Pp (PpfPp = —/). Whereas Pf and g ordinarily 
commute, when acting on the vortex core zero modes 
of the Af = 1 phase, these operators anti-commmute 
(PpgPpg = —1), since g flips the fermion parity of the 
states |±). 

This projective action of symmetry ensures that no 
symmetry preserving perturbation can gap out the vor¬ 
tex core zero-modes, and there is a two-fold symmetry 
protected degeneracy in every vortex core - indicating 
that the Af = 1 phase is an example of a projective f- 
SPT. 


2. Af = 2 , Fractional Symmetry Charge f-SPT 

We now review that the Af = 2 phase is a fractional 
symmetry charge f-SPT. Picturing the Af = 2 phase as 
a decoupled stack of chiral superconductors would sug¬ 
gest that a vortex contains four Majorana modes, includ¬ 
ing two Z 2 -neutral inodes: 7+, a =i,2 and two Z 2 -charged 
modes: 7_ ja =i,2- However, the 4-fold degeneracy associ¬ 
ated with these Majorana modes can be completely lifted 
without breaking symmetry, for example by adding the 
Z 2 symmetric perturbations: iy+qy+, 2 , and fy-yy-^. 

Despite having no symmetry protected degenerate 
modes, however, the vortex still has nontrivial symme¬ 
try properties. To see this, a useful mental device is to 
imagine promoting the Z^ fermion parity to a full 17(1) 
symmetry defined at the edge by: 

(lR/L,i\ u{ i) / cosd sin 6 \ h R / L ,i 
\"/R/L,2j y-sind COS ej \JR/L,2 


(and similarly in the bulk) and then later break this extra 
17( 1 ) symmetry back down to Zf. 

With the enhanced 17(1) symmetry, the Af = 2 phase is 
nothing but a bilayer quantum Hall (QH) state consisting 
of a v = 1 QH state of Z 2 -neutral particles and a v = — 1 
QH state of Z 2 -charged particles. This analogy, along 
with simple Laughlin flux-threading arguments, immedi¬ 
ately shows that the 7r-flux binds a half unit of Z 2 charge, 
i.e. transforms by a phase of ±* under the Z 2 symmetry 
(note the ± signs can be interchanged by locally binding a 
Z 2 -charged fermion to the 7 r-flux, but the integer charged 
fermions cannot screen the residual half-Z 2 -charge). 

This fractional charge is unchanged upon breaking 
the enhanced 17(1) symmetry back down to the origi¬ 
nal fermion parity, Zf. This implies that the vortex in 
the Af = 2 phase carries a fractional charge of the Z 2 
symmetry indicating that this phase is an example of a 
fractional charge f-SPT. 


3. Af = 4 , Bosonizable f-SPT 

Viewing the Af = 4: state as two copies of the Af = 
2 phase immediately shows that symmetry acts non- 
projectively on the vortex in the Af = 4 state, and that 
the vortex carries integer Z 2 -charge. In fact, we can al¬ 
ways consider the vortex to be Z 2 -neutral by combining 
it with Z 2 charged fermions to “screen” its charge (note 
that this binding of fermions does not change the bosonic 
statistics of the vortex, due to the mutual statistics be¬ 
tween fermions and the vortex). 

The trivial symmetry properties of the vortex imply 
that fermionic degrees of freedom do not play an impor¬ 
tant role in determining the topological properties of this 
phase - and that the fermion excitations can be gapped 
at the edge without breaking symmetry. However, the 
phase is not a trivial insulator. 

To see this, imagine inserting a symmetry flux (a tt- 
flux that is seen only by the Z 2 -charged fermions). As 
for the Af = 2 phase, it is again conceptually conve¬ 
nient to enhance the fermion parity to a full 17(1) ro¬ 
tation symmetry. Then, we may view the model as a 
stack of v = ±2 QH phase of Z 2 -charge ±1 particles. In 
particular, the symmetry 7r-flux binds integer Z 2 charge 
through the QH response. Such a 7 r-flux bound to unit 
charge object formally has semionic statistics (i.e. the 
edge operator corresponding to the symmetry flux will 
have semionic statistics, or equivalently, if we gauge the 
symmetry so that the symmetry flux becomes a gapped 
dynamical quasiparticle excitation, it will be a semion). 
As for the Af = 2 phase, this property of the symmetry 
flux is not altered when we dispense with the enhanced 
17(1) symmetry breaking it back down to the original Z^ 
fermion parity. 

This semionic statistics of symmetry-flux is the defin¬ 
ing property of the purely bosonic SPT with Z 2 symme¬ 
try introduced by Levin and Gu^. In fact, the Af = 4 
state is topologically identical to the Levin-Gu bosonic 
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SPT. Hence, we see that the Af = 4 phase is an example 
of a bosonizable f-SPT. 


4. Af = 8 , Trivial phase 


Using the now familiar tact of viewing the Af = 8 phase 
as two copies of J\f = 4 phases, we see that there are no 
nontrivial properties of either the vortex or Z 2 -symmetry 
flux in this phase. This implies that the A/" = 8 TSC is 
topologically trivial in the presence of interactions, and 
moreover that phases with Af differing by a multiple of 
eight are topologically equivalent. 

The remaining distinct phases with J\f = 3,5,6, 7 are 
combinations of other Af-phases described above, and 
hence are just combinations of projective and fractional 
charge f-SPT, and/or bosonic SPTs. 


Appendix D: Details of fermionic decorated domain 
wall construction 


2. Explicit form of “decoration” terms 


To explicitly write the decoration terms, it is conve¬ 
nient to represent the complex fermions by Majorana 
fields Xr,a = Cr t a + 4, a and Vr,a = ~i ( C r ,a ~ 4,a)) Such 
that c rjQ = (Xr,a + ir]r,a)/ 2. Eq. can be easily written 
as an explicit sum of local terms: 


ZJ _ (1 + T p)(l ~ r P')„._. 

^decorate — / J ^ ^\r,a^]r-\-d PP i ,a 


(PP'),r 


(1 — T I>)(1 + T P') 


2 

1 + T*pT Z p , . 


^Xr+dppi ,aVr,a T 
UG,a?7r,a (Dl) 


where P and P' are neighboring plaquettes, such that the 
line connecting P to P' bisects the bond connecting r and 
r + dpp '. An arbitrary convention must be chosen to fix 
the orientation of dpp': for example, we will take the 
Kitaev dpp> oriented clockwise around plaquettes with 
Tp, = — 1. The last term simply projects sites away from 
t— DWs into a trivial empty state. 


1. SPT Character of DDW model 


3. Odd fermion parity of Kitaev ring and related 
complications 


Suppose we can construct a Hamiltonian with symme¬ 
try group Z 2 x G, for which the DDW wave-function, 

I^ddw) = E c I^pt(C)) ® \t z (C )), is the unique 
ground-state. Here, C is a domain wall configuration cor- 

responding to r 2 wavefunction |r 2 (C)), and |V’g-spt 40) 
is a fermionic wavefunction with fermions along domain 
walls of C being in the ID f-SPT phase protected by sym¬ 
metry G, and fermions away from the domain walls be¬ 
ing in a trivial product state. Then, we can readily see 
that Tddw) represents a nontrivial SPT ground state. 
First, |'Pddw) preserves the G symmetry since r 2 trans¬ 
form trivially under G, and |V’G-Spt(C)) is manifestly G- 
invariant for all C. Second, since the r 2 variables are in a 
paramagnetic state |\&ddw) preserves the Z 2 symmetry. 

Lastly, I'Lddw) cannot be deformed into the trivial 

paramagnetic state |vf t riviai-PM) = Ec l^triv( C ))®l r *( C )) 
without breaking either Z 2 or G, where | v E , triv(C)) is a 
trivial product state of the fermionic degrees of free¬ 
dom. To see this note that, unless G is broken, upon 
deforming IV’G-Spt(C)) into |?4 ri v) one necessarily passes 
through a DW phase transition where low-energy exci¬ 
tations emerge on a DW of length Ldw with character¬ 
istic energy gap Ldw- Moreover, in order to preserve 
Z 2 the bulk wavefunction inevitably contains r 2 config¬ 
urations with DWs of arbitrary length, up to the system 
size L. Together this implies that deforming I^ddw) into 
l^triviai-PM) inevitably creates low-energy excitations ex¬ 
citations with energy ~ / that vanishes in the thermo¬ 
dynamic limit, i.e. implies that | v E , ddw) and |d'triviai-PM) 
are separated by a bulk phase transition. 


The TF terms written above face serious problems for 
Af 7 ^ 0,4. These issues stem from the fact that the 
ground-state of the Kitaev chain with periodic bound¬ 
ary conditions has odd fermion parity compared to the 
ground-state of a trivial chain of the same length. To see 
this, note that the total fermion parity of an L-site chain 
is given by: 

L L L 

P F = II(-1) C — = II *Xj\ra%a = i A.j,aVj+l,a 

3= 1 3=1 3= 1 

(D2) 

and recall that the trivial chain ground-state has 
iXj,aVj,a = IVj, whereas the topological one has 
i Xj ,aVj+i,a = IVj. Here site indices j are taken mod¬ 
ulo L. 

Hence, for Af = 1, the operators: 4/ c,p = \ p c(P))( P c\ 
appearing in Eq. [ 6 ] are fermionic whenever the parity of 
the number of domains domains differs between C(P) and 
C. These operators are not self-hermitian and square to 
zero: 4 /q P = {Fc\Fc(p))'&C,p = 0 , and anticommute on 
distant plaquettes: '&C,p'} = 0 (P/ P'). Hence, 

such terms are highly non-local and cannot appear in the 
Hamiltonian. This problem can be surmounted by imple¬ 
menting local constraints that prevent the global parity 
of domain walls from changing, however this necessarily 
spoils the exact solvability of the model. Therefore the 
DDW Af = 1 construction fails. 

For general Af, Z 2 domain parity changing terms take 
the form of Oc,at ~ flail ^c.p.a- For Af odd, O c ,u 
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are again fermionic and the construction fails just as for 
Af = 1. For even AT, on the other hand, all terms in Eq. [6] 
are bosonic. Despite this, the construction still fails for 

AT = 2,6. 

To see this, consider the Af = 2 case. To ensure that we 
obtain a nontrivial SPT phase the G symmetry must pre¬ 
vent terms like *xV,iXV ,2 and iVr,i r lr ,2 that can be added 
to adiabatically deform the Ising domain walls from dec¬ 
orated to trivial without closing a gap. However, the 
operators Oc ,2 ~ IIa=i 2 ^ e,p,a have the same form, and 
must therefore also be symmetry forbidden. Hence, the 
Af = 2 DDW construction inherently breaks the G sym¬ 
metry, and fails to produce a nontrivial SPT phase. Sim¬ 
ilar arguments rule out Af = 6. 

These problems are absent, however, when Af = 4. 
Here, it is possible to ensure that all trivializing bilinear 
terms (e.g. i\r,aXr,b ) are odd under symmetry, with¬ 
out ruling out the four-fermion operators GcyX neces_ 
sary to implement the DDW construction. Consequently, 
the DDW construction succeeds at providing a commut¬ 
ing projector Hamiltonian with nontrivial SPT character 
only for Af = 4. One suitable choice of a symmetry group 
G that can protect four copies of the Kitaev chain, which 
we will explore in the next sections, is to take G = Z|, 
and choose each of the four Kitaev chains to have a dif¬ 
ferent G-charge. 


4. Symmetric commuting projector character of 
Af = 4 model 

In this Appendix, we explicitly verify the commuting 
projector structure of the bosonizable Af = 4 TSC de¬ 
scribed in the main text. First, it is straightforward to 
verify that Eqs. [4] and [6] can by recast in the form of a 
commuting projector Hamiltonian with a unique ground 
state. The terms in Eq.[4]for fixed {P, P' , r}, and those in 
Eq.[6]for fixed P square to unity, and hence have eigenval¬ 
ues ±1. Then, for any such term, G, Ho,± = \ (1 ± O) 
are projection operators having eigenvalues {0,1} and 
satisfying ± = n o,±- Secondly, it is easy to verify by 
inspection, that all terms respect the required Zf sym¬ 
metry. 

It then remains to check that the projectors commute 
with one another. By inspection the decoration terms in 
Eq .[4]commute with each other. Also, by design, the dec¬ 
oration terms also commute with the terms in Eq.[6j since 
the transverse field terms preserve the low energy deco¬ 
rated DW subspace of ITdecorate> and vanish when act¬ 
ing on improperly decorated DW configurations. Lastly, 
starting with any spin configuration C and fermion state 
| Fc) and flipping two different plaquettes, P,P', using 
the terms in Eq. [6j produces the same spin configura¬ 
tion, and the fermion configurations that agree up to a 
potential overall phase: | Fc(PP')) = e*“| Fc(P'P))- Since 
the terms in Eq.[6]have the same action on all four flavors 
of fermions we may write e la = J | a=1 4 e * aa = ( eMl ) ■ 

To compute this phase, note that, since all terms in 


the Hamiltonian are real, the Hamiltonian satisfies an 
extra time-reversal-like anti-untiary symmetry that acts 
like complex conjugation T = K in the fock basis of the 
site fermions c r , a (he. Tc rj0 T _1 = c rja , Txr, a T^ x = 
Xr,a, Trir^T ^ 1 = —T) r ,a)- This, anti-unitary symmetry 
requires that a a G 0,7r, implying that e 4 * ai = 1, i.e. that 
the TF terms all commute. 

These considerations verify that the DDW construc¬ 
tion for the Af = 4 TSC protected by Zf symmetry has all 
the desired properties for a MB localizable parent Hamil¬ 
tonian. 


5. Equivalence of fermionic and bosonic phases 

for Af = 4 

In this Appendix, we show explicitly that the Af = 4 
TSC phase with Zij symmetry is equivalent, with inter¬ 
actions, to a bosonic SPT phase. 

Let us write the symmetry group as Z ^ 1 x Zf x Z % 
to label the different Z 2 subgroups. As described in the 
main text, a consistent, solvable DDW model of the Af = 
4 phase is obtained by decorating with four copies of 
Kitaev chains. A nontrivial SPT phase is obtained by 
taking each copy to have a distinct charge under G = 
Zf x Zj : {(+,+),(+, ^), (—, +)j (—j —)}• 

Without interactions, cutting open the decorating ID 
SPTs reveals four Majorana end states on each end - 
indicating 2 2 = 4-fold degeneracy. With interactions, 
we can gap out half of these degrees of freedom by 

adding the interaction term: —y(* 7 ++ 7 _|_ ){i'y _ |_7 _) 

that favors making the fermion parities P 2 = (?’ 7 ++ 7 h_) 

and P 2 = (*7 _ |_7 _) the same. There is still a two¬ 

fold degeneracy with this interaction term, spanned by 
|±) = |Pi = ±, P 2 = ±). Zf is diagonal in the |±) basis, 
and we may represent its generator as <m|±) = ±|±). On 
the other hand, we see that the Z^ generator changes the 
parity values: gsPi^gB = —Pi, 2 - Since gB interchanges 
| ±) it necessarily anti-connnutes with gc on the end-state 
manifold. This projective realization of the G symmetry 
protects the remaining two-fold degeneracy. 

The above considerations show, that the Z^ DW ends 

B C 

carry a “spin- 1 / 2 ” like degree of freedom on which Z 2 ’ 
act like 7r-spin rotations around orthogonal axes (e.g. are 
represented by the Pauli matrices: a x,v respectively). 
This is a discrete analog of Haldane’s famous gapped 
spin -1 chain, and is precisely the projective realization 
of Z| symmetry obtained by the purely bosonic ID SPT 
discussed in Refs. ITU and l28l 

Consequently, we see that, with interactions, this ID 
chain is topologically equivalent to the bosonic Z 2 x Z 2 
chain plus a trivial gapped fermion insulator. Conse¬ 
quently, the corresponding Z 2 x G DDW construction 
is topologically equivalent to the bosonic SPT state de¬ 
scribed in Ref. [28] plus a decoupled trivial fermion insu¬ 
lator. However, the fermionic and bosonic constructions 
are qualitatively distinct - since the fermionic DDW ap¬ 
proach taken in this paper produces a state with highly 
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entangled bosonic and fermionic degrees of freedom. We 
will see at the end of the following section that this quali¬ 
tative distinction makes the fermionic DDW state closer, 
in some sense, to the noninteracting topological super¬ 
conductor than to the purely bosonic one. 


6. Edge properties 

In this Appendix, we work out the edge properties of 
the MB-localizable Af = 4 TSC phase protected by Z| 
symmetry, and verify that they have the expected SPT 
properties. For clarity, let us distinctly label the three 
different Z 2 subgroups of the Z| symmetry as A, B, and 
C such that the Ising spins t z transform under Z^, and 
the fermions c r>a transform under Z^ x Z^\ 

To introduce an edge, consider the above Af = 4 model 
with spatially uniform couplings in a finite sample with 
open boundary along the x direction. In this section we 
identify possible edge phases of the the ground-state of 
the spatially uniform Af = 4 DDW model. These SPT 
edge state properties will be inherited by every state in 
the spectrum for the strongly random MBL version of 
this model model. 


a. Gapped symmetry broken edge states 

One can represent the edge degrees of freedom consists 
of complex fermions c rA on the edge sites, and Ising spins 
on T z r+1 on links between the sites on the edge. At low 
energies, only the MZM end states of r DWs are relevant. 
In terms of bulk fermions c r ^ a the DW zero modes are 
either \r,a or i] rta for t z _ 1 r T z r+1 = ±1 respectively. 

A trivially insulating edge can be obtained by plac¬ 
ing the boundary spins in an ordered phase with long- 
range correlations in rf r+[ . at the expense of break¬ 
ing the Z 2 symmetry at the boundary. Alternatively, 
we can produce gapped edge that preserves Z 2 but 
breaks time-reversal by choosing the edge spins to 
be eigenstates of r* r+1 , and gapping the low-energy 
fermionic zero modes with a Z^ symmetry breaking term: 
As Yja.b [iXr,aM a ' b Xr,b + M a:b r] r , b \ , where M ab is any 

suitable mass matrix. 

However, an edge that respects the full Z 2 symmetry 
must have Z^ domain walls proliferated, and the residual 
Z® x Zj' symmetry protects some degeneracies of the cor¬ 
responding DW MZMs. Hence, the proliferation of DWs 
containing protected zero modes, required to preserve Z^ 
symmetry, necessarily produces a gapless edge. 


b. Fermionic gapless symmetry preserving edge 

Insight into the gapless symmetry preserving edge ter¬ 
minations of the Af = 4 model can be obtained by consid¬ 
ering an antiferromagnetic arrangement of the Ising spins 


t z = (—1)"' where j labels sites along the edge. This ter¬ 
mination breaks Z^, but leaves a residual symmetry of 
7*2 = Zwhere T x generates a spatial translation by 
one lattice spacing along the edge (which we take to be 
along the x direction). 

In the absence of any edge couplings, there are four 
unpaired Majorana fermions 7 s ,s';j at each site along the 
edge for this AFM spin arrangement, where s, s' £ {±1} 
label the charges under Zf x Zff, such that the edge is 
highly degenerate. We may largely lift this degeneracy 
by adding the symmetry allowed couplings: 

_ iv ^ 

rJF.edge — ^ y 'Ys,s';j+l'ls,s , ;j 

j,s,s' 

= V ^ Sink 7 s , s ';fc 7 s , s ';-fc (D3) 

s,s';k> 0 

where 7 s , s ' ; fc = JT e 1 ^"/ s , s ';j- Linearizing the disper¬ 
sion near k = 0, tt edge mode in terms of right (R) and 
left (L) movers: 7 S s';R/L,j = e~ lk ^'y ss '-k that trans- 

/ 7T 

form oppositely under the modified Z^ symmetry. We 
can then imagine quantum disordering the AFM edge 
structure to restore the original Z ^ 1 symmetry without 
gapping out the fermionic edge modes. 

These arguments show that the gapless symmetry pre¬ 
serving edge of the Af = 4 model consists of four copies 
of nonchiral gapless Majorana fermions, each carrying a 
distinct combinations of Z^ x Z ^ 7 charges, and verifies 
that the solvable interacting lattice model based on the 
DDW construction reproduces the desired edge proper¬ 
ties of its free-fermion analog. 


c. Bosonic Gapless Edge 

As shown in Appendix |D 5| above, interactions can re¬ 
duce the ID fermionic SPT used to decorate DWs to a 
bosonic one while preserving the Z| symmetry. To see 
how this works in terms of the continuum description of 
the free-fermion edge with four non-chiral gapless Majo¬ 
rana modes, it is useful to use a bosonized description of 
the gapless fermionic edge. From the 8 chiral Majorana 
fermion fields: 7 s ,s',l//Jj define 4 chiral complex-fermion 
fields that can be bosonized: 

i’s^/R = 2 (7 s,+,L/R + ils,-,L/R ) = e~ l<t>L l R ' a (D4) 
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In terms of the boson fields 4> a< a the Zj symmetry acts 
as: 


4>r,s — 4>R,s 

<PL,s -► 9L,s + n 

7j B 

4>r/l,+ — <t>R/L,+ 

7j B 

&R/L- —^ <t>R/L- + 7T 
4>R/L,s —~~^R/L,a 

(D5) 

The edge Lagrangian density in the absence of inter¬ 
actions is then: 

£ed ge = ^ [d x <j) I K IJ d t (j)j + d x cj)iv LJ d x (j)j] (D6) 

where I,J £ {(s = +1, R), (-1, R), (+1, L), (-1, L)}, 
K IJ is a diagonal matrix with diagonal entries 
{+1, +1, —1, —1}, and v IJ is a non-universal velocity ma¬ 


trix that depends on the edge hopping amplitudes (and 
is also affected by interactions). 


By analogy to Appendix D 51 above, the edge inter¬ 
action term A J dx cos ( ^Z a =L/R s=±i 0a, s(x)J is allowed 
by symmetry. When relevant (in the renormalization 
group, RG, sense), this term gaps out all fermionic ex¬ 
citations at the edge. However, for weak interactions at 
the edge (small A and v IJ nearly diagonal) the fermion¬ 
gapping term is RG irrelevant - and the fermionic edge 
described in the previous section is stable to interactions. 


This highlights a qualitative distinction between the 
fermionic DDW construction here, and the bosonic one 
of Ref. [21J though the J\f = 4 topological superconductor 
produced by the fermionic DDW is topologically equiv¬ 
alent to a purely bosonic SPT accompanied by a spec¬ 
tator trivial insulator of fermions in the presence of in¬ 
teractions. Since the fermionic DDW approach produces 
a state with strongly entangled fermionic and bosonic 
degrees of freedom its symmetry preserving edge state 
is naturally closer to the free fermion edge with gapless 
Majorana modes. 
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